STM Imaging of Flux Line Arrangements in the Peak Effect Regime 
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We present the results of a study of vortex arrangements in the peak-effect regime of 2H-NbSe2 by 
scanning tunneling microscopy. By slowly increasing the temperature in a constant magnetic field, 
we observed a sharp transition from collective vortex motion to positional fluctuations of individual 
vortices at the temperature which coincides with the onset of the peak effect in ac-susceptibility. 
We conclude that the peak effect is a disorder driven transition, with the pinning energy winning 
from the elastic energy. 

PACS numbers: 74.60.Ge, 64.60.Cn 
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It is well known that the critical current density j c 
in weakly disordered type II superconductors shows a 
sudden increase when the applied magnetic field H ap- 
proaches the upper critical field H C 2 ■ An explanation 
for this intriguing phenomenon, known as the peak effect, 
was suggested by Pippard Q already in 1969. It is based 
on the competition between the collective work done by 
the disorder (the pinning centers) on the vortex lattice 
(VL) and the elastic energy stored in vortex lattice de- 
formations due to the pinning H . The work depends lin- 
early on (H C 2 — H) while the elastic energy essentially is 
proportional to the shear modulus cqq of the VL, which 
behaves as (i? C 2 — H) 2 . Therefore, at some field close 
to H C 2 , the pinning will exceed the elasticity and the VL 
will accommodate to the random pinning potential, lead- 
ing to the sudden increase in j c observed experimentally. 
The discovery of the high temperature superconductors 
made it relevant to also consider thermal fluctuations as 
a third energy scale, which lead to melting by suppression 
of C66- This would then yield a peak effect as a precursor 
of the melting transition j|, ||, || . 

Experimentally, it is hard to probe which mechanism 
is predominant, because one needs information about 
structural changes of the vortex lattice positional order 
as a function of time. Small angle neutron scattering 
(SANS) yields structural information along the field di- 
rection averaged over the sample volume, i.e., about the 
amount of entanglement of the VL. SANS has been suc- 
cessfully used in the peak-effect regime @, |[ [| but it 
is resolution limited in the transverse direction. Lorentz 
microscopy ]To|, scanning Hall probe (SHP) pj , scan- 
ning SQUIDJl|, || and, recently, magneto-optics ex- 
periments 1 14 yield positional information on the scale 
of individual vortices, but sense the magnetic field dis- 
tribution, and only work in the low flux density regime, 
usually far below H C 2, where vortices are well separated 
and the vortex-vortex interaction is very weak. In this 
Letter we show that information on the time-dependent 
positional order can be provided by scanning tunneling 
microscopy (STM), which is uniquely able to access the 
necessary length scales and flux densities. We present 
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FIG. 1: Phase diagram of NbSe2 obtained from susceptibil- 
ity measurements on a single crystal, showing Hci(T) (solid 
line / circles), the peak of critical current H P (T) (solid line / 
plusses) and the onset of the peak H v \ (T) (dashed line / trian- 
gles). The inset shows the peak effect in the bulk critical cur- 
rent density using the method due to Clem jl^, [l7| for a flat 
sample in perpendicular field vs temperature at /io-ff = 1.75 
T. 



STM data of vortex structures obtained on a pure single 
crystal of 2H-NbSe2 at temperatures around 4.5 K (far 
below T c = 7.1 K) in the regime of the peak effect. We 
observed a sharp transition from collective vortex motion 
to positional fluctuations of individual vortices at a tem- 
perature which coincides with the onset of the peak effect 
in ac-susceptibility measurements. The state above the 
onset temperature is characterized by the fact that the 
shear modulus of the VL does not play a role and shows 
the properties of a pinned liquid. 

Fig. ^ shows part of the H, T-phase diagram of one 
of our crystals as obtained from magnetization measure- 
ments with the field along the c-axis (perpendicular to 
the cleaved surface). Also shown (inset) is the critical 
current density j c as function of temperature T in a con- 
stant field of 1.75 T, obtained from the maximum of the 
out-of-phase component of the ac-susceptibility jl(| p7[ . 
With this method the same pinning conditions hold for 
the entire sample and geometrical effects are taken into 
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FIG. 2: Raw and processed data from one sequence of STM 
images acquired at 1.75 T for three different temperatures. 
The frame size is about 410x460 nm. Left (a-c) : original 
images (the dark elongated stain in the upper-right corner is 
caused by a surface defect); middle (d-f) : result of convolu- 
tion with the shape of "single vortex", right (g-i) : result of 
convolution with "perfect" lattice cell of 3x3 units. 



account. The peak effect sets in at T p \ ~ 4.5 K, it reaches 
a maximum at T p ~ 4.6 K, and superconductivity disap- 
pears at T C (H) ~ 4.7 K. It should be noticed that in 
dc-magnetization (SQUID) measurements no peak effect 
was observed due to fiux creep. This means that even 
in the peak-effect regime nonuniform flux distributions 
decay within a few minutes to a uniform flux density. 

The STM setup was the same as in a previous report 
ijhSf , with the sample immersed in liquid helium. After 
decreasing the magnetic field from above H C 2 at A.27K to 
a preset value (e.g. 1.75 T) in persistent mode (relative 
decay less than 5x 1CP 8 s -1 ), we waited about 10 minutes 
before starting the STM imaging. By closing the evapo- 
ration valve during imaging a slow temperature rise was 
established at a rate of 5xl0 -5 K/s to a value of 4.9 K. 
The temperature was determined from the pressure with 
a precision of 2 mK, but we found a small systematic dif- 
ference between the temperatures of the STM Tstm and 
the SQUID T S tm, T S q, T S tm -Tsq ~ 60 mK, due to 
the non-uniform temperature in the liquid. We corrected 
for this difference and estimate the residual uncertainty 
in the temperature determination to be about 20 mK. 
The bias voltage was set inside the superconducting gap 
(0.2-0.8 mV) and the tunnelling current was about 8- 
15 pA. The signal-to-noise (S/N) ratio goes down when 
approaching T C (H), e.g., the relative conductance modu- 



lation between vortex core and vortex cell edge decreases 
from 0.25 (far below T pl ) to 0.02 (at T p ). Therefore we 
had to use a low scanning speed of about 20-70 s per im- 
age to acquire images of size 512 x 512 pixels. The scan 
area was about 500x500 nm; at a field of 1.75 T we could 
probe the spatial arrangement of about 150 vortices. The 
images show a small constant spatial distortion due to a 
slight asymmetry in the piezo scanner. By heating up 
in constant field the vortex structure is very weakly per- 
turbed by the almost negligible change of the equilibrium 
magnetization. The effect of the tunneling current is also 
expected to be negligible because it is very small (< 15 
pA) and it is directed parallel to the vortices. Moreover, 
there was no influence on the vortex arrangement when 
we changed the tunnel current from 7 to 25 pA. In this 
respect our experiments differ essentially from dc and or 
ac-magnctization or transport measurements which do 
affect the vortex configurations. 

Results of experiments at 1.75 T are shown in Fig. 2a- 
c. We selected three typical temperatures: (a) 4.3 K 
(well below T pl ), (b) 4.44 K (just below T pl ) and (c) 
4.53 K (between T pl and T p ). The first image shows an 
almost perfect vortex arrangement without large visible 
distortions. With increasing temperature we see an in- 
creasing shear distortion of the VL and the appearance 
of VL defects. With a further increase of T the signal 
from the vortices gradually decreases to the noise level. 
Several image-processing procedures were used to extract 
the positional information. The S/N ratio could be con- 
siderably improved by convoluting the original data with 
the images of perfect vortex lattice cells acquired at low 
T, using cell sizes up to 10x10 vortices. In Figs, ^d-f the 
results are shown of convolution with 1 vortex ('lxl'), 
in Figs. |^g-i of convolution with a 3x3 cell, both with 
respect to the corresponding data in the Figs. ||a-c. Us- 
ing the lxl-convolutions, the displacements of individual 
vortices could be followed up to 4.56 K, see below. The 
low contrast in Figs. ||e,f is most likely due to noise, but 
may also denote local vortex displacements. 

With the knowledge of the vortex positions, further 
analysis can be performed using the vortex motion. In- 
specting consecutive frames in the temperature interval 
between 4.3 K and 4.6 K in a field of 1.75 T, we notice 
a sudden change in the vortex mobility at 4.48 K. Below 
4.48 K we find collective jumps of large portions of the 
VL; above 4.48 K, we only see oscillations of individual 
vortices around fixed positions. An example of a collec- 
tive jump is given in Fig. ||a, which displays the difference 
between two consecutive frames of raw data at 4.31 K. 
The scan direction is horizontal going from bottom to 
top. In the lower part of the figure the vortex lattice did 
not move substantially, and the difference gives an im- 
pression of the noise level. At the scan position denoted 
by the arrow the entire VL moved over a considerable 
fraction of the VL parameter do- The next difference 
picture (not shown), has a lattice below and noise above 




FIG. 3: (a) Difference between two consecutive images at 
4.31 K in a field of 1.75 T, to demonstrate coherent motion 
of vortices. When scanning the second image from bottom to 
top, the VL moved at the line indicated by the arrow, (b) Av- 
erage of the sum of 60 consecutive images taken with increas- 
ing temperature above T p i. The brightness is a measure for 
the probability of finding a vortex at a certain position. The 
numbers indicate individual vortices discussed in the text. 

the position of the arrow. We could see no net motion 
of the VL, it is jumping between metastable positions. 
Above 4.48 K, no collective jumps were observed; in the 
60 images in which T slowly rose to 4.9 K (total obser- 
vation time about 70 minutes) we did not detect any net 
motion of individual vortices in the image, although some 
vortices hopped back and forth between two metastable 
positions. A movie of the change in behavior, made of 
1-ccll filtered frames for noise reduction, can be viewed 
on our website |l9| . 

We can characterize this transition by introducing the 
following two order parameters, based on the displace- 
ment di^ n of vortex i going from frame n to frame n + 1, 
di.n = ^i.n+i^'Ti,^ (i) $n', in order to quantify the amount 
of motion of the vortices, we define 

^ := (^E d ln) 1/2 )/M, (1) 

i 

where the sum runs over the N displacement vectors, (ii) 
K n ; in order to quantify the amount of correlated motion 
of the vortices, we define 

ly 4n • d k n (2) 

where the sum runs over the displacements on neighbor- 
ing vortices i and k in a range 1.5 do- For both parame- 
ters, only contributions were retained for which 4« wa $ 
two or more pixels. 

Clearly, when vortices move coherently, both 5 n and 
K n will be large, while for uncorrelated, random motion, 
K n is small but 8 n remains finite. Both K n and S n are 
plotted versus T in Fig. ^a,b. For both, the behavior 
abruptly changes at T,sxm = 4.48 K, which is essentially 
T p i (see Fig. |l|). Below T pl the behavior of K n exhibits 
the collective vortex lattice jumps between metastable 
positions in the collective pinning potential. At T p i, the 
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FIG. 4: (a) Order parameter K„ for the correlated displace- 
ment of vortices from frame to frame as function of temper- 
ature at a magnetic field of 1.75 T. (b) Amplitude of vortex 
displacement S n averaged over the total number of vortices in 
a frame in the same sequences of images. Symbols (+) mark 
the value \/2 < u 2 >/ao (the value of <5„ for fluctuations of 
individual vortices) calculated from the average of 4 selected 
vortices (see Fig. pi and Table 1). 

picture suddenly changes to a situation in which each 
vortex independently fluctuates within its own pin po- 
tential. The sudden drop of K n to the noise level reflects 
this transition from collective to single vortex behavior. 
From the plot of 5 n versus T it follows that the collec- 
tive jumps span a distance of about 0.1 - 0.25 a , which 
changes at T v \ to individual jumps with a much smaller 
amplitude spread. The distance of collective jumps re- 
flects the typical length scale separating independent lo- 
cal minima in the pinning potential in the field regime 
close to H C 2 (at low fields where the vortex cores are well 
separated, this scale would be the Ginzburg-Landau co- 
herence length £). The jumps occur more frequently with 
increasing temperature. 

From the original STM pictures it can be also con- 
cluded that the jumps are directed along the principal 
axes of the VL. In the regime above T p i the mean- 
squared displacement gradually grows with temperature 
to a value <u 2 > 1 / 2 ^ 0.1 do- This number is close to the 
Lindemann criterion for melting, <u 2 > 1 / 2 = clclq with 
Cl ~ 0.15 — 0.25. However, we do not observe diffusive 
motion of vortices over more than eto. This is illustrated 
in Fig. 3b where 60 images in the temperature regime 
above T v \ are averaged. Since the noise level for all vor- 
tices is supposedly the same, the brightness of the spots 
represent the probability amplitude for finding a vortex 
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at a certain position. The brightest spots display vortices 
with an almost fixed position. One of them is marked #1; 
other vortices (marked #2 — 4) are less bright, indicating 
larger individual positional oscillations most likely caused 
by thermal fluctuations. 

TABLE I: Mean square displacement < u' 2 > x / 2 /a at three 
temperatures bei 
cated in Fig. 



Vortex number: 


#1 


#2 


#3 


#4 


4.51 K 


0.045 


0.095 


0.055 


0.077 


4.53 K 


0.055 


0.080 


0.13 


0.071 


4.55 K 


0.070 


0.10 


0.21 


0.15 



In Table 1 the values are given of <u 2 > 1 / 2 /ciq for 
these 4 vortices determined from 20 images around the 
denoted temperatures. Near the onset of the peak ef- 
fect (4.51 K) the amplitudes remain far below the Lin- 
demann criterion. Approaching T p , the amplitudes for 
the vortices #3 and #4 increase to values close to the 
Lindemann criterion. However, the estimated average of 
<it 2 > 1 / 2 /ao for all detectable vortices coincides with the 
data in Fig. Qb and remains well below the criterion for 
melting. Therefore, and because of the absence of large 
scale vortex diffusion within the 70 minutes observation 
time, we have shown that the change at T p \ signifies a 
disorder driven transition from collective pinning to sin- 
gle vortex pinning. The consequence is a sudden change 
in the temperature and/or field dependence of the critical 
current which causes the peak effect. 

Following this main result we add a few quantitative 
remarks about the observed disorder and the relation of 
our result to those in other papers. From the convoluted 
images we can extract some information about the size 
of the Larkin domains in the temperature range below 
the peak effect (although the limited number of vortices 
in the images does not allow much statistics). The posi- 
tional order in an area of 10x10 vortices disappears at a 
Tstm of about 4.35 K, while for an area of 6x6 vortices 
positional order is seen up to T pl . The computed corre- 
lation lengths from the Fourier transform of the images 
have similar values: (4 — 6)ao at the onset of the peak 
effect and ~ (1 — 2)ao above the onset temperature. It 
should be noted this the correlation length is about twice 
as large as the transverse size of the Larkin domain R c . 
At the onset of the peak effect the value of R c /a drops 
from about 3 to 1 . This result is in reasonable agreement 
with the values we obtained in a previous analysis of criti- 
cal current data on a similar crystal Jl7j . Based on SANS 
experiments on a Nb single crystal, Ling et al. || found 
evidence for supercooling and superheating of the vortex 
matter disorder around the peak effect line in the phase 
diagram and concluded that this is direct structural evi- 
dence for a first-order vortex solid-liquid transition. Very 



recent SANS work of Forgan et al. Q demonstrated that 
the vortex solid in very pure Nb is stable to very close to 
the H c2 (T) line. This shows that the history dependence 
observed by Ling et al. is disorder driven rather than of 
thermal origin and this is in agreement with the disorder 
driven transition seen in our STM experiment. 

In conclusion, the mechanism for the peak effect is 
the conquest of the pinning energy over the elastic en- 
ergy, just as was first suggested by Pippard [§), and 
later in more detail described by Larkin and Ovchinnikov 
f^of . The high temperature phase is characterized by 
the fact that the shear energy has become irrelevant. In 
this phase, the vortices are not free to move around as 
one would expect to occur in a liquid, but they remain 
trapped in their own pin potential which may differ for 
each vortex. In terms of the collective pinning theory this 
would be the amorphous limit, but it would be equally 
correct to say that we are dealing with a pinned liquid. 
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